arXiv:1503.03937v3 [hep-ph] 19 Nov 2015 


Prepared for submission to JHEP 


FERMILAB-PUB-15-079-E 


Dark Matter Directionality Revisited with a High 
Pressure Xenon Gas Detector 


Gopolang Mohlabeng," Kyoungchul Kong,^ Jin Li,^ Adam Para,^ Jonghee Yoo^ 

^Department of Physies and Astronomy, University of Kansas, Lawrenee, KS 6604-5, USA. 

^ Center for Underground Physies, Institute for Basie Seienee (IBS), Daejon 305-811, Korea 

^Fermi National Aeeelerator Laboratory, Batavia, IL 60510, USA. 

E-mail: gopolang.mohlabeng@ku.edu, kckong@ku.edu, jinlee@ibs.re.kr, 
para@fnal.gov, yoo@fnal.gov 

Abstract: An observation of the anisotropy of dark matter interactions in a direction-sensitive 
detector would provide decisive evidence for the discovery of galactic dark matter. Directional infor¬ 
mation would also provide a crucial input to understanding its distribution in the local Universe. Most 
of the existing directional dark matter detectors utilize particle tracking methods in a low-pressure 
gas time projection chamber. These low pressure detectors require excessively large volumes in order 
to be competitive in the search for physics beyond the current limit. In order to avoid these volume 
limitations, we consider a novel proposal, which exploits a columnar recombination effect in a high- 
pressure gas time projection chamber. The ratio of scintillation to ionization signals observed in the 
detector carries the angular information of the particle interactions. 

In this paper, we investigate the sensitivity of a future directional detector focused on the pro¬ 
posed high-pressure Xenon gas time projection chamber. We study the prospect of detecting an 
anisotropy in the dark matter velocity distribution. We find that tens of events are needed to exclude 
an isotropic distribution of dark matter interactions at 95% confidence level in the most optimistic 
case with head-to-tail information. However, one needs at least 10-20 times more events without head- 
to-tail information for light dark matter below ^50 GeV. For an intermediate mass range, we find it 
challenging to observe an anisotropy of the dark matter distribution. Our results also show that the 
directional information signihcantly improves precision measurements of dark matter mass and the 
elastic scattering cross section for a heavy dark matter. 
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1 Introduction 

The Standard Model (SM) of elementary particle physics has been astonishingly successful in ex¬ 
plaining much of the presently available experimental data. However, it still leaves open a number of 
outstanding fundamental questions whose answers are expected to emerge in a more general theoretical 
framework. One of the major motivations for pursuing new physics beyond the SM is the ‘dark matter 
puzzle’, which finds no explanation within the Standard Model. From the accumulated experimental 
data, we now know that ordinary matter comprises only about 4.9% of the Universe. The remaining 
95.1% is divided between a mysterious form of matter called ‘dark matter’ (26.8%) and an even more 
perplexing entity called ‘dark energy’ (68.3%) [1]. 

Naturally, discovering dark matter (DM) and measuring its properties has become central to 
the fields of particle physics, astrophysics and cosmology. The diversity of possible dark matter can¬ 
didates requires a well-balanced program based on direct detection experiments, indirect detection 
experiments, collider experiments and astrophysical probes sensitive to the non-gravitational interac¬ 
tions of dark matter. Vast experimental and technological progress in the coming decade will put the 
most promising ideas to the test [2]. 

In the standard scenario of a WIMP (weakly interacting massive particle), direct detection ex¬ 
periments record the nuclear recoil energy spectra produced when a dark matter particle scatters off a 
target nucleus. The expected nuclear recoil energy falls exponentially and such events, with an energy 
of typically not more than a few tens of keV, lie well within the range of abundant backgrounds due to 
radioactivity and other cosmogenic backgrounds. Despite these challenges, experimental limits on the 
interaction cross-section versus WIMP mass have been steadily improved. Yet, there exists no widely 
accepted evidence of their presence on Earth. Firm evidence of directionality relative to a WIMP wind 
would be the most robust signature of the WIMP nature of dark matter, and is an essential step to a 
claim of discovery. If any of the direct detection experiments observe evidence for nuclear recoils that 
cannot be explained with known processes, then the search for directionality in such nuclear recoils 
will be of foremost interest. Given that dark matter signatures exhibit an exponentially falling energy 
spectrum, and an annual modulation of interaction rates which can be easily mimicked by any activity 
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of seasonal variations, or cosmogenic backgrounds, a discovery claim of dark matter would have to be 
followed by other exceptional evidence [3-6]. 

A powerful signature of dark matter would be the detection of a significant spatial anisotropy 
in the angular distribution of such nuclear recoils consistent with the standard model of a non-co- 
rotating WIMP halo. Earth’s position in the galactic arm provides a boost of ^230 km/s, comparable 
to the quasi-virial velocity ^220 km/s of gravitationally captured WIMPs. Dark matter interactions 
would produce a large forward-backward asymmetry in the angular distribution of nuclear recoils. An 
importance of the directional information for an incontestable claim of discovery and its role in the 
additional rejection of terrestrial backgrounds is widely appreciated [7-17]. 

Currently most attempts at directional detection have focused on low-pressure gas time projection 
chambers (TPCs), in order to provide the 3-D track reconstruction and energy resolution needed to 
identify low energy nuclear recoils [18-29]. In all low-pressure TPC detectors, strong tension exists 
between the desire to use a very low gas density so that nuclear recoil tracks are long enough to be 
imaged with adequate clarity, and the desire to increase the gas density so that greater sensitivity 
can be realized. The Diffusion of the ionization image during drift, limits the drift length and the 
avalanche amplification noise and/or photon detection quantum efficiency degrades the quality of the 
track information. In most cases, the total mass per detector in these approaches is less than a kg 
due to the limited scalability of the low pressure detector. Thus progressing to ton-scale masses would 
imply a very large and impractical number of separate devices. 

Recently, a novel approach has been proposed to confront the challenges of the directional sensi¬ 
tivity for nuclear recoils with active masses approaching the ton-scale. The detector concept is based 
on a high-pressure Xenon gas TPC with an electroluminescent gain stage which utilizes the ‘columnar 
recombination’ (CR) process, leading to a potential directional sensitivity of nuclear recoils [30]. If 
this conceptual idea and the related detector technology can be demonstrated, it would revolutionize 
dark matter experiments [31]. Unlike low pressure gas tracking detectors, a ton-scale high-pressure gas 
dark matter detector would be more practical. 

In this paper we examine the sensitivity of the proposed high-pressure Xenon TPC directional 
dark matter detector. We pay special attention to the capability of measuring head-to-tail information 
and distinguishing the incoming and outgoing direction of the recoiled nucleus. We study the prospect 
of detecting an anisotropy in the dark matter velocity distribution. We begin our discussion with a 
short review on directional dark matter detection in Section 2. We devote Section 3 to a detailed 
analysis involving a high-pressure Xenon gas detector. We discuss more on columnar recombination, 
set up our numerical study (Section 3.1) and examine how much improvement can be made on the 
measurements of mass and cross-section (Section 3.2). We then investigate the angular distributions 
and anisotropy of the dark matter distribution in Section 3.3 and Section 3.4. Section 4 is reserved 
for discussion. 

2 A Brief Review on Dark Matter Directionality 

The motion of the solar system relative to the Galactic WIMP halo provides a distinctive signal for 
WIMP detection. This circular orbit of our solar system around the galactic center results in a very 
strong forward-backward asymmetry in the angular distribution of nuclear recoils produced in WIMP 
events. The differential nuclear recoil rates as a function of both recoil energy and recoil angle have 
been extensively studied in literature. Recoil rates including the angular distribution of events were 
first discussed in Ref. [7] and then further developed in Refs. [3, 4, 32]. In this study we adopt the 
formalisms used in Refs. [3, 4, 33] and only provide a short review in our paper. 
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Figure 1. Geometry of a WIMP scattering off a target nucleus in the detector. The WIMP is incident 
at an angle (a, P) relative to the z axis. The nucleus recoils in the (^, (j)) direction. 


Let us consider a WIMP particle of mass incident at velocity v = v{sm a cos ;dx+sin a sin /3y-\- 

cosaz) in the detector, as illustrated in Figure 1. After interaction with a WIMP, the target nucleus 
recoils with some velocity u = u{sm 0 cos (jyic + sin6>sin0y + cos6>z) and momentum q at a direction 
{O^cj)). The rate at which this nucleus recoils per unit recoil energy per unit recoil angle is given as 
follows: 
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where po is the dark matter halo density in our local part of the galaxy, aw n is the WIMP-nucleus elas¬ 
tic scattering cross-section, and F{q) is the nuclear elastic scattering form factor. Assuming that the 
nucleus can be approximated to be a sphere with uniform density, the form factor is the Fourier trans¬ 
form of the nuclear density. This gives us the Helm form factor F{q) = ^i^'^AQrn)-qrnCos(qrn)] ^ 

where is an effective nuclear radius and q = \j2M^Fji is the recoil momentum of the nucleus. Fji 
is the recoil energy of the nucleus and r = 4 :MnM^/{Mn + ^ kinematic factor [32]. Angle 0 

is the recoil angle which determines the direction between the recoiling nucleus and the initial WIMP 
trajectory and /(v) represents the velocity distribution of WIMPs in the galactic halo. We call Eq. 
2.1 the double differential recoil rate. Mn = 0.932A GeV is the target mass, with A the atomic mass 
number of the target atom in atomic mass units (AMU). The factor 0.932 is the value of AMU in GeV 
and pn = M^M ^/+ M^) is the reduced mass of the WIMP-Nucleus system. 

The double differential recoil rate can be defined in a simpler mathematical form by adopting the 
following mathematical convention [3]: 


/K,q) = 




■q - Vq)f{-V)(f 


( 2 . 2 ) 


where Vq is the minimum velocity a WIMP must have to impart a recoil momentum q to the nucleus, 
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or equivalently to deposit an energy Er = q the recoil momentum direction and v is the 
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velocity of a WIMP particle in the halo. Eq. (2.2) is the definition of a three-dimensional Radon 
transformation, which represents the velocity distribution for a stationary detector in the galactic 
frame. For an observer moving with velocity Viab in the galactic frame, this is the velocity of the 
observer in the galactic frame. It is related to the velocity of the WIMP in the lab frame, viab and in 
the galactic frame Vgai by a Galilean transformation viab = Vgai — Viab- The properties of the Radon 
transformation for a pure translation Viab [3] imply, 

flahiVq, q) = /galK + Viab-q, 4) • (2-3) 


For our study we assume a Maxwell-Boltzmann velocity distribution, truncated at the escape 
velocity Vesc of the WIMPs ; 


fM{v) 





(2.4) 


for V <Vesc and /m('^) = 0 otherwise, k^sc is a normalization factor which is obtained by integrating 
the velocity distribution in the galactic frame from 0 to Vesc- For the velocity distribution Eq. (2.4), 
the radon transform becomes. 
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Finally for a detector on earth moving through the galaxy with velocity ve in the direction of Cygnus 
X-2, Viab = Ve and q ■ Viab = —ve cos9. We can combine Eqs. (2.1) and (2.5) to obtain 
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where Nq is the Avogadro’s number and Vq = v^in = \/Er /with Eq = the most 

probable kinetic energy of the WIMPs. We choose the most probable WIMP velocity vq = 2‘^^km/s 
and the escape velocity of the WIMPs from the galactic halo, Vesc = 600/cm/s [32]. is calculated 
in the appendix of Ref. [32] and includes the velocity of the Earth with respect to the Sun, the 
proper motion of the Sun and the velocity of the solar system with respect to the galactic center. The 
WIMP-Nucleus cross section is defined as cfwn = E {A — Z)fn)^ where A = ^/jcrwp^ 

and similarly for /^, with fp and fn the WIMP-proton and WIMP-neutron couplings respectively. 

2 

In the case where fp ^ fn (which we assume), we obtain cfwn = ^cfqA^ with ctq = (Twn = crwp? 
the WIMP-nucleon cross-section and /ip ~ /i^ [32]. As the truncated Maxwell-Boltzmann velocity 
distribution exhibits a rotational symmetry (along 0), Eq. (2.6) is only dependent on the polar angle, 
0 . 


3 Dark Matter Directionality with a High-Pressure Xenon Gas Detector 

In this paper we focus on two aspects of directional detection, “parameter estimation” and “measure¬ 
ment of anisotropy”, with emphasis on a high-pressure Xenon gas detector. As such a detector does 
not currently exist and only a concept is discussed [30] , it is uncertain what features would be appro¬ 
priate to consider. Therefore we assume certain important detector parameters for our study. We first 
introduce four different types of detectors (Section 3.1) for discussion and compare their performance 
in the parameter estimation of the WIMP mass and WIMP-nucleon cross section (Section 3.2). Fur¬ 
ther information on the angular distributions and an anisotropy in the dark matter flow are presented 
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in Section 3.3 for those different detectors. Although they are generic dark matter detectors without 
details of a particular detector concept, results are still relevant to grasp potential performance of a 
high-pressure Xenon gas detector and how it would compare as to non-directional detectors. Finally 
in Section 3.4 we study anisotropy of the WIMP velocity distribution using a high pressure gas TPC. 
Throughout our studies, we include detector resolution with relevant energy threshold cuts. 

3.1 Columnar Recombination and the Numerical Set Up 

The importance of directional information in dark matter experiments has been recognized for a long 
time [6, 19-21, 34-37]. The short range of the low energy nuclear recoils is an obvious experimental 
challenge. Therefore a low pressure gas TPC is a natural experimental choice to extend the observable 
length of the recoil trajectory up to macroscopic dimensions thus enabling the determination of the 
spatial direction of the recoil track. Unfortunately about 1/10 bar of low pressure gas limits the 
practically attainable mass of the detector. Hence the technique is currently only applicable in case 
of relatively large interaction cross sections. 

Recently a conceptual dark matter detector that utilizes columnar recombination (CR) [38] has 
been proposed by D. Nygren [30] as a possible technique for the determination of the recoil direction 
in massive detectors, up to several tons. The detector exploits CR within an ensemble of ions and 
electrons generated by the nuclear recoil. The CR process occurs when the detection of a highly ionizing 
track and an externally applied electric field coincide, such that the external field drives the ionization 
electrons to drift in the vicinity of the ion column. Electrons traveling at distances close to the ions 
and below the Onsager radius undergo electron-ion recombination with an emission of characteristic 
photons. Conversely, recombination is much less likely if the particle track and the electric field are 
perpendicular. When the particle track and E-field coincide maximum CR is expected as opposed 
to when they are parallel to each other. The amount of CR can give us an estimate of the relative 
angle between the track and the E-field, thus a measurement of the nuclear recoil angle. Hence, the 
directional information of nuclear recoils might be obtainable in a high-pressure gas detector. 

This preliminary detector concept explores the possibility of utilizing a special Penning mixture 
in the Xenon gas which will convert the energy harbored in primary excitations to ionization. In 
addition to the uniform drift electric field and a charge collection plane it would be equipped with 
internal reflectors and photodetectors allowing for a highly efficient collection of light from the entire 
volume of the detector. After the initial interaction, the electrons drift along the field direction towards 
the collection plane. The electrons may undergo recombination, with the emission of characteristic 
photons. The number of the emitted photons and thus the size of the light signal S', will depend on the 
angle between the recoil track and the drift field direction, therefore it is suggested that the division 
of a total signal, S + J, into its components S and ionization I should depend on the angle Ol between 
the recoil track direction and the direction of the electric field E in the TPC. 

Practical implementation of such a concept awaits experimental demonstration [31, 39]. In par¬ 
ticular the head-to-tail capabilities of the detector are of great importance. The direction of recoil is 
determined through cos Ol = fig^) therefore its values are limited to be positive. If the columnar 
recombination is forward-backward symmetric it will allow for the determination of |cos^l|, otherwise 
the mapping of S/{S + I) onto |cos^l| may be multi-valued, but there will be some region around 
cosOl = 1 (most likely) or cos Ol = —1 with characteristically higher values of S/(S + /). In the latter 
case one will be able to classify all events into two angular bins: larger or smaller than some cos6>o, 
where a possible value of Oq must be established experimentally. 

We compare the physics potential of various classes of detectors with different capabilities of the 
directional measurement for studies of the dark matter interaction with a cross section of awn = 
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5 X 10“^^ pb = 5 X cm^. We also assume that the forthcoming generation of experiments will 
focus on the demonstration of the galactic origin of the observed signal (directionality) and on the 
determination of the properties of the dark matter (interaction cross section and mass). Our analysis 
is restricted to dark matter signal events assuming zero-backgrounds, which is a good estimate for a 
4 keV threshold cut 

To be specific, we consider the following progression of possible detectors: 


• a detector with no directional capabilities 

• an ‘ideal’ detector capable of measuring the recoil angle in the range — 1 < cosO < 1 , with some 
characteristic resolution (we refer to this case as ‘head-to-tail’) 

• a ‘symmetric columnar recombination’ detector, thus capable of the determination of |cos^l|, 
with some characteristic resolution (we refer to this case as ‘no head-to-tail’, it is also known as 
a ‘folded’ directional rate [41]. See also Ref. [42] for related studies.) 

• an ‘asymmetric columnar recombination’ detector, thus capable of classifying the events in two 
angular bins 


In addition we compare the capabilities of directional detectors constructed on a movable system 
that maintains the orientation of the detector’s electric field in the galactic frame - thus rotating in 
the Earth coordinates. We use the direction of the Earth’s motion as our reference direction, and we 
define the direction of dark matter flow as our forward direction, which is opposite to Cygnus. We 
call this the ‘parallel’ case, when E is aligned with our forward direction and ‘perpendicular’, when 
E is perpendicular [43]. We define the corresponding angle between the electric held and the recoil 
direction as Ol = 0\\ and Ol = respectively. This set up conveniently identifies Ol = 0\\ as the 
recoil angle 0 in Eq. (2.6), for the parallel case. 

Eigure 2 illustrates the double differential distributions for the first two types of detector concepts 
(head-to-tail or no head-to-tail) for both a parallel and perpendicular electric held. We define them 
as “ (a). dgJlLeni “ (b), (c) and “ (d), respectively, where 

subscripts || and T denote the direction of the drift electric held with respect to the WIMP direction. 
Along each curve, the same number of events are expected. Eigure 2(a) is the most ideal case with a 


full coverage of the recoil angle. By our set up. 


dN 


= n and also 


dN 


, which is 


dcos^ii dcosO ^|cos0||| d|cos0|‘ 

the ‘folded’ directional recoil rate, where | cos^| does not distinguish the beginning of the recoil track 
from its end (lack of head-tail discrimination) [24, 41]. A detector that is fixed on Earth may weaken 
the DM directionality and we have investigated this effect in our simulation by orienting the electric 
held at a fixed angle a with respect to the incoming WIMP direction as shown in Eigure 3. The 
(a = 0 case corresponds to a movable detector that we have described and the detector that is fixed 
on Earth would include a combination of different a angles, washing out the angular information. As 
the movable system provides the best sensitivity, we will consider this case throughout the paper. 

Note that our study point 5 x 10“^^ pb for a light dark matter particle falls within the over¬ 
whelming neutrino backgrounds in direct detection experiments as described in [40]. The effects of 
neutrino backgrounds on directional detection have been partially studied in Ref. [39]. 

To determine the dark matter mass, cross section, and anisotropy, we perform simulations for these 
types of detectors. We assume an energy threshold of 4keVnr (unless noted differently). Gaussian 


^Very little neutrino background is expected for recoil energies above 4 keV [40], e.g., 
events for a Xenon detector with 10 ton-year. 


about 0.5 neutrino 
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(b) Parallel, No Head-Tail 



(c) Perpendicular, Head-Tail 



(d) Perpendicular, No Head-Tail 
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Figure 2. Contours of constant number of events in a two dimensional space for a 10 ton-year Xenon 
detector. Events are normalized to a case with an 80 GeV dark matter particle and a WIMP-nucleon 
cross-section of 5 x 10“^^ pb. We consider a detector with (a) head-to-tail capability and a parallel 
electric field ^ (b) no head-to-tail capability with a parallel electric field (c) 

head-to-tail capability with a perpendicular electric field ( cieS^s e± ) head-to-tail with a 

perpendicular electric field 


smearing is applied for both energy and angle as follows: 


F{E,e) = I F{E',e') I" 




■ 2ai 
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(3.1) 













COS0|| |cos0|| 

Figure 3. Angular distribution of nuclear recoil events with (a) and without (b) head-to-tail capability, 
rotated at different angles with respect to the incoming WIMP trajectory for a 80 GeV WIMP, a = 0° 
and a = 90° have been discussed earlier. 


where F{E,0) is event rate function (Eq. 2.6), aE = A\/E is the energy resolution and ao is a 
constant angular resolution. We have assumed A = 1 and ag = 30° in our numerical study, unless 
noted otherwise. They are rather conservative choice compared to those reported in literature (see Ref. 
[6] for details.). In the case of low energy recoils we would have to worry about negative energies in the 
above Kernel, but we found that a threshold cut at 4 keVnr is large enough to avoid such events. The 
angular smearing was carried out in 6>'-space using the Kernel in Eq. 3.1. Eor a given number of events 
at an angle O-hin (0 < 6> < tt), the smearing Kernel is applied to a large array of linear angles-bins of 
0'. The events that fall below 0 and above tt respectively are then folded back on the main range of 
the distribution. This is done to preserve the angular range of the original 0 distribution and in this 
way the total number of events is conserved as required. We choose a cross-section of 5 x 10“^^ pb 
for simulation purposes (unless noted otherwise), which roughly gives 103 (143) events after (before) 
the 4 keV threshold cut in a Xenon detector for 10 ton-year, assuming zero background. 


3.2 Estimation of WIMP Mass and WIMP-Nucleon Cross-Section 


To illustrate the capability of measurement of the WIMP mass and cross section for the detectors 
discussed in the previous section, we calculated expected event rates using Eq. (2.6) for four different 
study points = 20 GeV, = 60 GeV, = 80 GeV and = 100 GeV for a fixed input 
cross section of awno = 5x10“^^ pb, as shown in Eigure 4. The event rates are normalized to a 10 ton- 
year exposure of the Xenon detector. The physics information (dark matter mass and interaction cross 
section) is determined using a binned likelihood analysis assuming a Poisson probability distribution 
of our signal events: 


^bin 

^=n 


i=l 




exp 


-Nl 


(3.2) 


where is the expected number of (template) events, Ne = NE{M^^awn) and Nq is the observed 
number of (signal) events in each bin. 














Figure 4. Illustration of measurements of mass and cross section of galactic dark matter at 90% 
confidence level (CL) with crwno = 5 x 10“^^ pb for (a) = 15 GeV, (b) = 60 GeV, (c) 

= 80 GeV, and (d) = 100 GeV. The four different contours represent different types of 

detectors assumed in the likelihood analysis. For ‘Energy0Angle’ (black, solid), we use both recoil 
energy and angular information obtained from the theoretical distribution. We integrate over the angle 
and annual modulation to obtain the recoil energy information only for ‘Energy only’ (red, dashed). 
A detector without head-to-tail information is shown in the green-dot-dashed contours. Einally for ‘2 
Angle Bins’ (blue, dotted) we use 2 bins in the angular distribution. The input point for our simulation 
is (M^O: cTwno) and is represented by a dot inside the ellipses. All events are normalized to 10 ton-year 
exposure for a Xenon detector, including detector resolution effects and 4 keV threshold cut. 


Eour different detectors are examined: (i) a conventional non directional detector where only the 
recoil energy of the events is measured (denoted as ‘Energy only’) shown as red-dashed ellipses in 
Eigure 4, (ii) a detector which has the ability to measure the energy, angle and annual modulation 
signal of every event (denoted as ‘Energy0Angle’) shown as the black-solid ellipses, (iii) a detector 
without head-to-tail information (shown in green-dot-dashed ellipses), and (iv) a detector in which 
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we do not have the ability to measure the angle of each event, but we can determine the number of 
events within a certain angular ‘cone’, i.e. we split the angular distribution of events in two bins, 
and use both bins in the likelihood analysis, but since we cannot determine the precise angle of each 
and every event, we only know that an event fell in this angular space of a certain size (with ±30° 
opening angle for this study) or outside (2 bin angle system). The ‘2 Angle Bins’ case in blue-dotted 
ellipses represent the performance of this type of detector. For representation purposes, we do not 
mention annual modulation in the figures, since the results do not change much whether the annual 
modulation effect is included or not. We carry out the binned likelihood analysis and obtain a region 
of the parameter space that is consistent with the input point at 90% C.L, shown as four ellipses for 
each case in Figure 4. The minimum of the log-likelihood is marked for each case and they should 
coincide with the input study point in the absence of any statistical fluctuations. Although finite 
statistics would shift the best fit point off from the original input and may alter the shape of contours 
slightly, our study indicates what improvement is expected in the best case scenario. 

The difference in the orientation of the ellipses for the 20 GeV case and the rest is easily understood 
from the interplay between a threshold cut and the shape of the differential energy rate. For a light 
DM of mass (20 GeV in this case), the differential distribution of recoil energy is very steep and a 
majority of the events are cut away with a threshold cut, which implies that one needs a higher cross 
section to fit the data with On the other hand, the fitting procedure requires a smaller 

cross section, as the energy distribution is less steep for This is shown in Figure 4(a). 

However this is no longer true, if the input mass is large as illustrated in Figure 4(b)-(d). 

We also notice that for a light DM, directional information does not play an important role in 
measurements of parameters. However for a heavy dark matter (heavier than 100 GeV), the full 
directionality is crucial in precision measurements. Precision can be substantially improved for the 
intermediate mass range below 100 GeV even without head-tail information. We illustrate this in 
Figure 5 where the relative WIMP mass uncertainty (SM^/M^) is shown as a function of WIMP mass 
for different classes of detectors. 

We have also studied an impact of both angular resolution and energy resolution to see the effect 
of finite resolution on the dark matter parameter determination and we find that the 90% contours 
do not change significantly. Also in the case of a detector with the 2 bin angular system described 
above, we used a benchmark angular opening of 30° in Figure 4. We then tested for different sizes of 
angular area, i.e. 60° and 90°, but we found no large difference in neither the angular distributions nor 
the 90% GL contours. Our results imply that directionality with full angular coverage improves the 
measurement of masses and cross section significantly especially for a heavy dark matter. A detector 
without head-to-tail information or one with limited 2 angular bins provides a marginal improvement 
in the accuracy of the parameter determination. 

3.3 Angular Distributions and Anisotropy of Dark Matter Flow 

A non-trivial angular dependence of nuclear recoils produced in dark matter interactions arises due 
to an asymmetric velocity distribution of dark matter in the lab-frame. As mentioned in the previous 
section, we use the truncated Maxwell-Boltzmann distribution in Eq. (2.4) as our reference. Therefore 
in this set up, the asymmetry shows up entirely due to the motion of the detector as shown in Eq. 
(2.6), i.e., is isotropic (flat) for ve = 0. 

To maximize the observed anisotropy we constantly adjust the orientation of our detector with 
respect to the Cygnus direction. As an exercise we have studied two detector configurations, one with 
the drift field parallel to the dark matter direction and one with the drift field perpendicular to it. In 
Eigure 6, we show angular distributions of two types of detectors based on the capability of head-to-tail 
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Figure 5. Relative WIMP mass uncertainty as a function of WIMP mass. We compare three detectors, 
the conventional (non-directional) detector (in red), one with full coverage of directionality (in black), 
and the other without head-to-tail information (in green). The case with 2 angle bins lies between the 
conventional case and the no head-to-tail case and is not shown here. A 15° (30°) angular resolution 
is assumed for the dotted curve (solid, both green and black). 


discrimination for several choices of dark matter mass. The two plots in the top panel correspond to 
, in (a) and i in (b), which is n \ , also known as the ‘folded’ differential in Ref. [411. 

A precise comparison would be somewhat difficult since we assume a different set up and different 
materials than those used in Ref. [41] (CS2 and CF4). However we are able to reproduce a (roughly) 
consistent result and especially the shape of our folded distribution resembles that in Ref. [41]. In our 
set-up, this folded distribution is obtained when the drift electric field is parallel to the initial WIMP 
direction based on details of the columnar recombination effect. , and , are shown in (c) 

and (d), respectively. 

As shown in Figure 6(a), dark matter scatters predominantly in the forward direction and details 
of the shapes of the angular distributions are dependent on the mass of dark matter and the imposed 
threshold cut. Unfortunately a lack of head-tail discrimination puts a severe hurdle on the measure¬ 
ment of the anisotropy of the dark matter flow (see (b) and (d)). The situation may be (slightly) 
improved for a light dark matter (20 GeV and 60 GeV) at the cost of signal statistics by imposing a 
higher threshold cut as shown in Figure 7, while there is no change for a relatively heavy dark matter 
(100 GeV). This is due to a correlation between the recoil energy and the scattering angle. Angular 
distributions observed in a detector with the field perpendicular to the average WIMP direction also 
show a similar behavior but the sensitivity of the measurement would be greatly reduced. 

3.4 Anisotropy with a Columnar Recombination Detector 

Our discussion in the previous section is somewhat generic in a sense that results do not particularly 
utilize the effect of CR. In a real experiment, the GR detector would not measure recoil energy and 
angle directly, but would rather count the number of electrons and photons released from Ionization 
(/) and Scintillation (S) processes respectively, with some detector resolution and efficiency. The recoil 
energy and the recoil angle are obtained as a function of the two variables I and S. The efficiency 
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Figure 6. Comparison between a detector with capabilities of measuring head-to-tail (left) and one 
which is incapable of such a distinction (right), for a 4 keV energy threshold cut, assuming the Drift 
Electric field is pointing parallel (top) or perpendicular (bottom) to the initial WIMP trajectory. 


of the measurement of / and S will also depend on the orientation of the drift E with respect to 
the nuclear recoil trajectory. If the field is parallel to the recoil trajectory, one expects a higher rate 
of ionization and scintillation, and the opposite effect when the field is perpendicular to the recoil 
trajectory. After these are measured one can convert these observables to recoil energy and recoil 
angle. We assume the following relation. 

S = F{Er) e Er cos^ Or , (3.3) 

I = F{ER)eEnsmHL, (3.4) 

where E{Er) is the number of observed photo-electrons per keV which takes into account the quenching 
factor, and e is the detection efficiency of photons. For the E{Er), we have adapted experimental 
values for absolute SI (prompt scintillation) yields for electron recoils in Xenon as in Ref. [44]. 
Absolute yield (in photons/keV) is given as a function of the incident gamma energy compared with 
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Figure 7. Angular distributions for a drift electric field oriented parallel (top) and perpendicular 
(bottom) to the initial WIMP trajectory for a detector with no head-to-tail capabilities. These are 
illustrated for 4 keV, 10 keV and 20 keV energy threshold cuts for 20 GeV, 60 GeV and 100 GeV 
WIMP masses. 


their Monte Carlo output taking into account the recombination probability. We have taken the best 
reproduction based on their model and scaled it down with our choice of e. Smearing S and / with 
Poisson statistics is performed before converting back to energy and angle. 

We simulate the amount of scintillation light and ionization yield that would be obtained for 
several WIMP masses with a cross-section of 5 x 10“^^ pb in 10 ton-year exposure of a Xenon detector. 
Converting (S', I) to {Eji, \ cos^^l), we obtain the results shown as the solid histogram in Figure 8 for 
two different e’s. The dotted histograms are illustrated in Figure 7, generated from an MG based on 
the theoretical expectations, assuming Gaussian smearing with 4 keV energy threshold. It is crucial 
to achieve a high photon detection efficiency as shown in Figure 8(b), which shows a close match 
with results (dotted histogram in Figure 8) in Figure 7. However at a dark matter mass of 80 GeV, 
the I cos ^ 11 1 distribution becomes flat even with a 10 ton-year exposure. The peak at cos^y = 0 is a 
result of the very low light signal expected in these conditions and Poisson fluctuations shifting events 
towards S' = 0, hence cos^y = 0. It gets reduced for a better efficiency and/or a higher threshold cut 
as shown in Figures 8 and 9. With a higher energy threshold cut, hence higher photo-statistics, the 
correct form of the angular distribution is recovered. 

We perform a likelihood ratio test to compute the required number of events to rule out the 
hypothesis of an isotropic velocity distribution. The results are shown at 95% CL in Figure 10. The 
test has been done for three cases: (a) head-to-tail case, (b) no head-to-tail capability, (c) the case 
with two angular bins. Solid (dashed) curve represents the exclusion limit (at 95% CL) that rules 
out hypothesis of a flat angular distribution (isotropic dark matter distribution) with 4 keV (10 keV) 
energy threshold. In other words, one can expect to observe anisotropy of dark matter distribution in 
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Figure 8. Angular distribution, for a realistic CR detector for two different photo detection 

efficiencies e. The solid histograms are angular distributions obtained from the ionization and scin¬ 
tillation light assuming a detector resolution with Poisson distribution using, cos Ol = \/S/{S 1). 

The dashed histograms are events generated from a Monte-Carlo simulation based on the theoretical 
distribution, assuming a detector resolution in Eq 3.1. 





Figure 9. Effect of a threshold cut in angular distributions. 4 keV (a), 10 keV (b) and 20 keV (c) 
energy threshold cuts are applied in each panel with a 20% detection efficiency for 20 GeV (in black), 
60 GeV (in red), 80 GeV (in blue) and 500 GeV (in green). 


the parameter space above the exclusion curve. The numbers in red give the required number of events 
needed at a certain WIMP mass to rule out an isotropic dark matter flow for 10 ton-year high-pressure 
Xenon detector exposure. Results are presented for an electric field parallel to the direction of the 
WIMP particles. We have imposed 30° angular resolution. The results are obtained with a theoretical 
double differential distribution, consistent with the dotted histogram in Eigure 8. The shaded area 
represents the current exclusion limit by LUX [45] after an appropriate rescaling of 0.0275 ton-year 
exposure to our case (0.0275 ton-year = 118 kg x 85 days). 
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Figure 10. Test of the anisotropy for (a) head-to-tail capability, (b) no head-to-tail case, and (c) the 
case with two angular bins. The solid (dashed) curves represent the boundary of exclusion of a flat 
angular distribution at 95% CL with 4 keV (10 keV) energy threshold. The numbers (in red) provide 
the estimated number of events required to exclude a flat angular distribution, i.e., for directionality 
at a certain WIMP mass. An electric field parallel to the recoil trajectory is considered for 10 ton-year 
high-pressure Xenon detector exposure. The shaded area represents current exclusion limit by LUX 

[45]. 


4 Discussion and Outlook 

An observation of the anisotropy of dark matter interactions would provide decisive evidence for the 
discovery of dark matter [13-16]. We investigated the feasibility of a high-pressure Xenon TPC dark 
matter detector which is sensitive to the angles of recoil produced in the interaction of dark matter 
particles with nuclei in the detector. The angular information helps precision measurements in the 
parameter space of the cross-section vs WIMP mass. Our study shows that full angular coverage and 
directionality could significantly improve the precision of the determination of the dark matter mass 
and/or the interaction cross section, especially for a heavy dark matter. The improvement is marginal 
for a detector without head-to-tail information or for a detector with the 2 angle bins. We find also 
that angular resolution does not make much difference in the improvement of the DM signal. 

The angular information of the recoil trajectory can be used to establish the anisotropy of the 
observed signal. A Xenon detector with a 10 ton-year exposure with head-tail capabilities could 
demonstrate the anisotropy of an observed signal if the corresponding interaction cross section is of 
the order of tens of 10“^^ pb. The sensitivity of a detector without head-tail capabilities would be 
reduced by at least one order of magnitude. It is interesting to notice that precision measurement is 
sensitive to a heavy dark matter while the anisotropy probe is more sensitive to a light dark matter 
for this type of a high pressure Xenon detector. 
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